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Abstract: We study the equilibrium measure for a logarithmic potential in the 
presence of an external field K(0 +^ PiO^ where t is a parameter, K(0 is a smooth 
function and p{^) a monic polynomial. When p{^) is of an odd degree, the equilibrium 
measure is shown to be supported on a single interval as \t\ is sufficiently large. 
When p(^) is of an even degree, the equilibrium measure is supported on two disjoint 
intervals as t is negatively large; it is supported on a single interval for convex p{x) 
as t is positively large and is likely to be supported on multiple disjoint intervals for 
non-convex p{x). 

The support of the equilibrium measure shrinks to isolated points as \t\ +oo in 
all the cases that we consider. For sufficiently large \t\, each topological component 
of the support contains a local minimizing point of the external field V^{^) + t p{C,)', 
a "potential well" phenomenon. 



§ 1 Introduction 

In this paper, we study the following minimization problem with constraints 



Minimize 

{V>0,/Vd?=l} 



1 

'2^ 
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The external field V{^) is a C°° function that satisfies 

Vie) 

lim , = +00 . (1.2) 

i^±^ log{l + 

Under this condition, the existence and uniqueness of the minimizer for (jl.lj) has been 
established [20] • The measure ip{^)dC,, where is the minimizer of is called 

the equilibrium measure under the external field V{^). 

Equilibrium measures find applications in many branches of mathematical sci- 
ences. It is used to describe the partition function of the Hermitian one-matrix model 
in random matrix theory of statistical physics |21IIH1- It is also intrinsically connected 
to the free energy of the Yang-Mills theory ^Hj- Finally, it plays an important role 
in orthogonal polynomials and approximation theory fK)\ . 

Although its importance to physics and approximation theory (0 1201) ^^^^ 
known, the minimization problem (ll.lj) is not well understood. The minimizer is 
explicitly known only for a few cases where the external fields V{C,) are the simplest 
polynomials [71 120] . It is therefore desirable to study the minimization problem for 
much more general C°° external fields. 

Our method to solve the minimization problem (jl.lll is the same as one [l2j we 
used to solve a similar minimization problem for the zero dispersion limit of the KdV 
equation 

Ut + QuUx + ^^Uxxx = with u{x, 0; e) = uo{x) , 

where the initial data uo{x) is a bounded decreasing function. The weak limit of the 
KdV solution u{x, t; e) as e — is determined by a minimization problem [T3lll6| l^ 



Minimize 

{ip>0, i/jSLi} 



1 

'2^ 



log 



HOHv)d^dr]+ / V{^,x,t)ij{Od^ 
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where the space-time dependence is through the function V 

V{^,x,t)=x^-Aet-9iO (1.4) 

and 6{^) is encoded with the initial information Uq{x). We found that the problem 
fjl.3|) is intrinsically connected to the Euler-Poisson-Darboux equations We uti- 
lized the solution of the equations to study the minimizer and hence the weak limit 
of the KdV solution. 

The minimizer of (ll.lj) has been known to have a compact support if the external 
field V{^) satisfies condition ()1.2|1 . The support is usually a union of a finite or infinite 
number of disjoint closed intervals. In this paper, we will study how the number of 
gaps in the support of the minimizer varies with respect to the external field. We 
would like to know whether there is any universality on the number of gaps. 

We shall consider parameter-dependent external field V{^) = K(0 + where 
t is a parameter, K(0 is a smooth function of ^ and p(^) a monic polynomial. This 
is motivated by the space-time dependence of V of ()1.4|) in the KdV minimization 
problem p.Hj) . We shall show that when p{^) is of an odd degree, the minimizer of 
p. HI has no gap in its compact support as \t\ is sufficiently large. When p(^) is of an 
even degree, the minimizer has a single gap in its support as t is negatively large; it 
has no gap for convex p{^) as t is positively large and is likely to have multi-gaps for 
non-convex p{^). We note that the same minimizer may have an arbitrary number of 
gaps when |t| is not large. A similar result has also been discovered in the case of the 
KdV weak limit [22: . Namely, the KdV weak limit is of zero or single phase for all x 
when t is sufficiently large, i.e., the minimizer of (jl.lll has either zero or one gap in 
its support. 

The support of the minimizer shrinks to isolated points as \t\ —>■ +00 in all the 
cases that we consider. Each topological component of the support contains a local 
minimizing point of the external field V{C,)] a "potential well" phenomenon. 

The organization of the paper is as follows. 

In Section 2, we will use function theoretical methods to solve the minimization 
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problem. We will formulate the minimizer in terms of solutions of the Euler-Poisson- 
Darboux equations. 

In Section 3, we will use the Euler-Poisson-Darboux solutions to study the behav- 
ior of the minimizer when the parameter in the external field is sufficiently large. 

§ 2 Solution of the Minimization Problem 

In this section, we will use the method that we have developed for the KdV zero 
dispersion limit iI2 to solve the minimization problem 
Introducing a linear operator 

L^(0 = -/ \og\^ - ^Im^^)d^^ , (2.1) 

71" J-00 

we rewrite the quadratic functional of as < L^, ^ > + < ^ >. Here 

< > is the standard inner product. The Euler-Lagrange equations take the form 

L^(0 - V^(0 = I where ^/^ > , (2.2) 
LipiO - ViO < I where = , (2.3) 

where / is the Lagrange multiplier. It can be shown 6] that ip is the minimizer iff is a 
nonnegative function that satisfies variational conditions ()2.2II2.3|1 and the constraint 



+00 

= 1 . (2.4) 

oo 



We make the ansatz that the support of ip consists of a finite union of disjoint 
intervals. One denotes J = ip > and writes 



/ = uf;i(M2fc,M2fc-i) , (2.5) 

where ^2^+2 < • ■ • < M2 < ^^i- Hence, the support is the closure of /. 
We now consider a slightly stronger version of ()2.2|) and ()2.3|) . 

^7/^(0- no = / on/, (2.6) 

Lm-V{i) < I onM\/, (2.7) 



where / denotes the closure of I. Since ■i/' > 0, we must also have 

^ > on / , (2.8) 
^ = off J . (2.9) 

Our strategy to construct the minimizer is to first find the solution ip of equations 
f|2.4|) . ()2.6|) and ()2.9|) and then impose inequalities ()2.7|) and ()2.8p on ip. Since it is a 
non- negative function and since it satisfies ()2.2|1 . ()2.3|1 and ()2.4|1 . will then be the 
minimizer. 

To solve ()2.4|) . ()2.6p and ()2.9p . we use the fact that the operator L of ()2.1|) is 
connected to the Hilbert transform H on the real line. 

[^Hij{T)dT, (2.10) 
Jo 



where 

V'(/^)c^/^ 

This makes equations (j2.6|) and (j2.9|) amenable to the Riemann-Hilbert technique in 
function theory. 



Hi^iO = -P.V. f 

J-c 



Differentiating ()2.6|) with respect to ^ and using ()2.1Up . one obtains 

= V'{i) on I . (2.11) 
To recover equation ()2.6|) from ()2.11|) by integration, one must have for /c = 1, 2, ■ ■ ■ , 

/ [H^ljiO-V'md^ = ■ (2.12) 

Recalling the relation of the Hilbert transform to analytic function, one can write 
for real ^ 

where g~^{C,) is the boundary value on the real axis of a function 

1 /-^"^ ^(/x)d/x 

g{z) 
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which is analytic in the upper half complex plane. In view of the constraint ()2.4|) . we 
expand the Cauchy integral to obtain Q{z) = —-7=^ + 0{l/z^) for large \z\. 
Conditions ()2.12|1 now take a new form 

[Img+{^)-V'{OH = Q (2.13) 

for = 1,2, ■ ■ ■ 

Equations (|2.9|) and (|2.1H) then become a Riemann-Hilbert problem in function 
theory 

Img+iO = V'iO on I , (2.14) 
ReG+iO = off/, (2.15) 

where Q{z) is analytic in the upper half complex plane. It follows from the Plemelj 
formula that 



TT J J R{fl, U){fl — z) 



where u denotes (mi, U2, ■ ■ ■ , U2g+2) and i?(^, u) = a/(^ - Ui)(^ - U2) ■ ■ ■ - M25+2), 
is a solution to this Riemann-Hilbert problem. Here R{^, u) is set to be positive for 
> Ml. It defines a Riemann surface with branch cuts along the set / of (j2.5j) . 
To derive the equations governing the endpoints Mi, M2, ■ ■ ■, U2g+2, one usually 
imposes condition ()2.13|) and the asymptotics Q{z) = — + OiXj z^) for large \z\ 
The former gives g equations and the latter results in (7 + 2 moment conditions. 
Consequently, one obtains exactly 2(^ + 2 algebraic equations on ui, M2, ■ ■ ■ , M2g+2- 

In this paper, we take a slightly different approach. We observe that the Riemann- 
Hilbert problem 1)2.14112.13)) has many other solutions. Indeed, it is obvious that 

Giz) = " ^ ^ (2.16) 

i?(z, u) 

where Q{z) is an arbitrary polynomial with real coefficients, is also a solution. We 
choose polynomial Q{z) such that 

1. it is a polynomial of degree 2g + 1. 



2. g{z) = -^=^ + 0(1/^2) large \z\ 

3. conditions (12.131) are satisfied. 



Hence, constraint ()2.4|) and conditions ()2.6|) and ()2.9p are built in the construction 

ofg. 

It is quite obvious that such a polynomial Q is unique. 

We first analyze the boundary value of g{z) on the real axis. Necessarily, the 
Cauchy integral in the numerator of ()2.16|) will become a singular integral. Our key 
observation is that the latter is, more or less, the solution of a boundary value problem 
for the Euler-Poisson-Darboux equations. 

Proposition 2.1. The boundary value of Q{z) at ^ Ui, i = 1,2, ■■■ ,2g + 2 

g+^^. ^ -2V^R\^, u)^g{^, ui, U2,--- , U2g+i) + v^Q(0 ^ yZi^'(^) (2 17) 

R{^, u) 

where (^g{S,,ui, ■ ■ ■ ,^23+2) satisfies the Euler-Poisson-Darboux equations 



2(^-Ui)^-^ = ^-2^, (2.19) 

We will omit the proof here, since it is similar to the proof of an analogous result 
on the KdV zero dispersion limit |12j . 

A formula for Q can be derived from the loop conditions ()2.13|1 and asymptotics 
Q{z) = — + 0(1/2;^) for large \z\. To achieve this, we introduce a sequence of 
polynomials 

P,,„(e, u) = + ag^^e^^-' + ■■■ + (2.21) 
whose coefficients are uniquely determined by 

^^ = r-^ + 0(^) for large lei, (2.22) 



and 



d^ = k = l,2,---,g . (2.23) 
Proposition 2.2. 

2g+2 2g+2 ^ 
i=l l=l,l=/=i 

^ CgPg,g(^,U) . 

Here ^E'g(^, Ui, U2, ■ ■ ■ , ^29+2) satisfies the Euler-Poisson-Darhoux equations h2.1^) and 
\2.1f^) with the diagonal boundary value 

1 d^^^Viu) , 



r/ie coefficients 



g-k 

Cfc = 2A; ^ Ti{u)qg^k+i{u) k = 1,2, - ■ ■ ,g, 

1=0 



where Ti{u) 's come from the expansion 



u) = /i^+Mro(n) + liM + IlM + . . . ] . (2.25) 
/i /i^ 



The function Qg^k is 



= I ■ (2-26) 

27rv — 1 Ji P{t^, u) 



The proof is also similar to that of the KdV case |12j . 

The function qg^kiu) also satisfies equations of Euler-Poisson-Darboux type [T^ . 
Finally, we postulate the boundedness of the minimizer, which is equal to the real 
part of Q^{S,)- The numerator in ()2.17p must then vanish at the end points of set / 

P(u„m) = z = 1,2,---, 2(7 + 2, (2.27) 

where 

P{i,u) =2R\i,u)^g{i,u) - Q{i,u) . (2.28) 



Since it is of degree 2(yf + 1 in ^ and has 2g + 2 zeros because of ()2.27|) . the polynomial 
Q{^,u) must be identically zero. The real part of of ()2.17|) . on the solution u of 
equations ()2.27j) . becomes 



= -2ReW-lR{i,u)}^g{^,u) . (2.29) 

Hence, V(0 of satisfies conditions (jZll), and (ESI)- 

Inequality ()2.8|) is equivalent to 

Re{y/^R{^, u)}^g{^, u) <0 for ^ on I , (2.30) 

and ()2.7p is equivalent to 



/ -R(/i, u)^g{fi, u)dfi > for 'U2fe+i < ^ < U2k and 1 < A; < 25f + 2 , 

f R{n,u)^g{n,u)dn > for ^ > ui , (2-31) 

J ui 

R{fj.,u)^g{^,u)d^ < for ^ < 'U2c,+2 • 



U2g + 2 



We summarize the above results in the following theorem 

Theorem 2.3. // {ui,U2,-- - ,U2g+2) satisfies equations \2.21\) and if inequalities 
\2. yUlXK'^} are satisfied, then ipiO liKEW ^■^ non-negative function that satis- 
fies variational conditions \2.i^) and \2.'J\) and constraint \2.4^ ; so ipiO '<^ust he the 
minimizer of the minimization problem M.l]) . 



§ 3 Large Parameter Results 

In this section, we shall consider a one-parameter family of external field 

v{i) = v,{0 + tp{0 ■ (3.1) 

Here, t is a parameter, p{S,) is a monic polynomial of degree n and K(0 is a 
C°°(— oo, +oo) function. We are interested in the behavior of the equilibrium measure 
when |t| is large. 
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Our strategy for studying the large parameter behavior of the minimizer of (jl.ll) 
is as follows. For simplicity, we will assume K(0 have a power function growth 
at ^ = — oo or ^ = +00. This allows us to use the scaling technique to study the 
equilibrium measure when the parameter t is sufficiently large. 

§ 3.1 The degree of polynomial p(^) is odd 

We shall show that the equilibrium measure is supported on a single interval 
[m2,'Ui] when \t\ is large. This corresponds to ()2.5p for g = 0. 

We will use Theorem 2.3 to construct the minimizer of Hence, we need to 

solve equations ()2.27|1 for (7 = and verify inequalities ()2. 30112. 3l|) . 

We now consider the algebraic equations ()2.27j) for g = 0. The function P(^, ui, U2) 
of (IT^ is 

2(^ - ui){i - M2)$o(^, Ml, U2) + ^2)^0(^1, Ml, U2) + Ui)^o(m2, Ml, U2) 



TT 



Equations (|2.27|) exactly become 



(mi -M2)^o(mi,Mi,M2) = 0, (3.2) 

(m2 - Mi)^'o(m2,Mi,M2) = 0, (3.3) 

TT 



which are equivalent to 



(mi - M2)^^^o(mi,Mi,M2) - y ;^ = 0, (3.4) 

^'o(ui,Ml,M2) + ^o(m2,Mi,M2) = 0. (3.5) 

In the derivation of ()3.4j) . we have used an identity 

d 

^o(mi, Ml, M2) - ^'o(m2. Ml, U2) = 2(mi - M2)7^^'o(mi, Ml, U2) ■ (3.6) 

C'M2 



This is easily verified by calculating 
d_ 



d 

2(e - M2) — ^'o(e. Ml, U2) - ^o(e. Ml, U2) 
M2 



^2 ^ Q 

2(e - M2)7r^*0 + 27^^0 - T^^O • 
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The right hand side vanishes because satisfies the Euler-Poisson-Darboux equation 
(I2.19|) . The function in the parenthesis is then independent of ^; so we obtain 

d 

OU2 

Letting ^ = ui yields the identity ()3.6|) . 

We will rewrite equation ()3.5|) in another useful form. Its left hand side is a 
function of Ui and U2. Denote this function by H{ui, U2). Since ui, U2) satisfies 
equations ()2. 18112. we derive an equation for H{ui,U2) 

d^H _ dH dH 

duidu2 dui du2 

The boundary condition ()2.24j) for \E'o('C, ^i, ^2) implies H{u,u) = V'{u). We then 
use formula ()B.3|1 to obtain 

1 n v'{ij,)dfi 



H{ui,U2) 



U2 - U2) 



IT 

Hence, equation ()3.5p is equivalent to 

r. (3,7) 

>y{ui - - U2) 

We shall study the case when t is negatively large in detail. The other case when 
t is positively large can be handled in the same way. 

We now make assumptions on Since V{^) = K(0 + ^ P(0 satisfies the 

growth condition p.2|) . 14(0 must grow faster than p{^) as ^ ^ +00 because of 
t < 0. In contrast, we will make a very mild assumption on K(0 ^ ^ ~1 

V40 = C^^^'^ + h^O fore»l, (3.8) 

K'(0 < fore«-l. (3.9) 

Here (7+ is a positive constant and is a positive constant that is bigger than n 
and 2. The number n is the degree of the polynomial p{^)- Function /i+(C) has an 
order less than n. More generally, we assume 

hm ^ = . (3.10) 
10 



This immediately implies 

ltoMi)^„_ Hm^ = 0, lim^ = 0. (3.11) 

We will solve equations ()3.4|) and ()3.5p for large ui and M2. This is motivated by 
the fact that V{^) has a critical point between U2 and mi in view of equation ()3.7j) 
and that the minimizing point of V{^) = K(0 + P{0 moves to +00 as t ^ —00. 

We will first consider the case p{^) = 

We will split \E'o of ()3.4II3.5|) into simpler terms. In view of its boundary data ()2.24|) 
for g = depends linearly on V. The decomposition of ^ = t^'^ +C+C,^'^+^^ + h^{^) 
allows us to write 

^0 = + + ^h+ii^ . (3.12) 

In view of the integral formula ()B.4|) . \I'^n and \E'^a/+ are homogeneous functions of 
Ml, U2) of orders n — 1 and M+ — 1, respectively. 

M_,_ 

Substituting ()3.12|) into equation ()3.4|) and then dividing it by \t\'^'-'^'+~"\ we use 
the homogeneity of \E'gn and \£'^m+ to obtain 



5 ...... ^ d 



2 



Id 1 1 
+ M+^2 —^h+{o{uuUi,U2)}2 M+ = • (3-13) 



Here ?7i and U2 are 



1^1 M^—n M^—n 



Similarly, substituting ()3.12j) into equation ()3.5|1 and then dividing it by \t\'^'+~" , 
we obtain 

-'^^n{Ui,Ui,U2) + C+'^^M+{Ui,UuU2) + i^-i ^fe+(g)(m,m,M2) (3.15) 

-^e(f^2,t/i,f/2) +C+^^M+(f/2,t/i,f/2) + i^,, ^^+(0(^2,^1,^2) = . 
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It follows from the behavior ()3. 10113. iTj) of that the three terms involving 

()3. 13113. 1'Sj) . together with their first derivatives with respect to Ui and U2, 
decay to zero as t —00 if Ui and U2 of ()3.14j) are kept bounded. 

We now denote 1/t by T and solve equations ()3.13|) and ()3.15|) for Ui and U2 as 
functions of T in the neighborhood of T = 0. 

First, equations ()3.13|) and ()3.15|) at T = become 



(f/l-t/2)-y/-^Vl/^n(f/i,f/i,f/2)+C+^Vl/^M+(?7i,f/i,f/2) = 0, (3.16) 

Uu U2) + (f/i, f/i, f/2) 

-^e(f^2,t/i,f/2)+C+*^M+(f/2,f/l,f/2) = 0. (3.17) 

These two equations have a solution Ui = U2 = U* , where U* is determined by 

where we have used the boundary condition ()2.24|) for \E'^n and \E'^a/+. We hence 
obtain 

Tl 1 

solution of (13. 16113.1711 . 

Second, we calculate the Jacobian of equations ()3.13j) and ()3.15j) at Ui = U2 = U* 
and T = 0. Denote the left hand side of ()3.13|1 by Fi{Ui,U2,T) and that of ()3.15|1 by 
F2{Ui,U2,T). Since \E'g satisfies the Euler-Poisson-Darboux equations ()2.18|) . ()2.19p 
and ()2.24|) . we obtain 



2 -e + c+e^+]i5=t7*>o 
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Hence, the Jacobian of ()3.13|) and ()3.15p is nonzero. Equations fj3.13|) and ()3.15|) then 
give Ui and f/2 as functions of T near T = 0. 

Therefore, equations ()3. 2113 .31) have a solution ui(t), U2{t) for negatively large t. 
The solution has the following asymptotics 



uiit) ^ U*\t\^'+~'' , U2(t) ^ f/*|t|"+-" for t < -1 . 

Moreover, it follows from ()3.13p and ()3.14j) that the length of the interval [ui{t), ui{t)\, 
the support of the minimizer, shrinks to zero as t —00. 

To make sure that of ()2.29j) is the minimizer, we also need to verify that 
inequalities ()2. 30112. 31)) are satisfied. 

We split $0 in the fashion of ()3.12p 

Using the scaling ()3.14|) . we write ^Q{^,Ui{t),U2{t)) as 

1 

\t\ (S, f/i, f/2) + C+$^M+ (S, f/i, f/s) + wi, U2)] , (3.18) 

M_|_— n 

1 

where S = We use (j3.11|) to deduct that if S is kept bounded, say, 

f/*/3 < S < 3U*, the last term in the parenthesis goes to zero uniformly as t ^ —00. 
The sum in the parenthesis then has the limit 

u*, u*) + c+<i>^M+ (s, u*, u*) 



2(S - t/*) 

2(S - U*) 



= ^ pi i>0 for^<H<3f/*. 

2(H - [/*) 3 - - 

Here we have used formula ()B.5|) for and in the first equality and U* = 

[n/C_|_M+]^/'^*''^+~") in the last two equalities. In view of (|3.18|) . we have the inequality 
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^o{^,ui(t),U2{t)) > for ^U2{t) < C, < 2ui{t) when t is negatively large. This proves 
the inequahty ()2.30p . 

To prove inequality ()2.31|1 . we first show that (^Q{^,Ui{t)^U2{t)) > uniformly for 
< — 1 and ^ > 2ui{t) when t is negatively large. 
We shall prove the uniform positivity for ,^ < — 1 first. We use formula ()B.6|) to 
write ^o{^,ui(t),U2(t)) as 

1 



V(^-^l)(^-^2) 



1 r V(e-wi)(e-w2) v'{fi)dfi 



TT 



U2 



^ - v/(ui - - U2) 

no 



_ 1 r W(e-^i)(e-^2) _ v^'(/i)rf/i 



(3.19) 



where in the last equality we have used equation ()3.7p . 

We now show that the integral of ()3.19|) tends to zero uniformly for ,^ < as 
t ^ — 00. Because of the inequality 

\- T 1| < < for ^ < < ^2 < Ml , 

4 - /i M2 - 4 M2 

the integral is bounded by 



U2 Ju2 a/(Mi - U2) 

for ^ < 0. This, in view of the scaling ()3.14|) . is further bounded by a constant times 

Plit) - U2{t)l 



Ml —n 



U2{t) 

M+-2 

where we have used V'{^) = 0{{ui{t) - U2{t))t^'+~") for U2{t) < ji < ui(t). This 
asymptotics can be verified by observing that V'{^) has a zero between U2{t) and 
uiif) because of ()3.7|1 and by expanding V'{n) around the zero. Since Ui{t) — U2{t) = 
0{\t\ 2{M_|_-n)^ because of ()3.13|) . the integral of ()3.19|) tends to zero uniformly for 
^ < as t — > —00. 

14 



The first term V'{C,) = tn^"'~^ + Vl{^) in the parenthesis of ()3.19|) is bounded from 
above by a negative constant for ^ < — 1 when t is negatively large. This immediately 
follows from condition 

Therefore, ^Q{C,,Ui{t),U2{t)) > uniformly for ,^ < — 1 when t is negatively large. 

In the same way, we can prove that $o(^) Ui{t),U2{t)) > uniformly for C, > 2ui(t) 
when t is negatively large. This together with $o(^, 'Ui(t), 'U2(i)) > for ^U2{t) < ^ < 
2ui{t) proves the first half of (IZ!r[|l . 

It remains to prove the rest of ()2.3H) . i.e., 

U2{t) 

^/ifi-ui{t)){fi-U2it)) ui{t),U2it))dfi < (3.20) 

for ^ < U2(t). Since $o(/i, 'Ui(t), 'U2(^)) > uniformly for ^ < — 1 and ^ > U2{t)/2 
when t is negatively large, it suffices to prove the inequality for — 1 < < U2{t)/2. 

Using formula ()B.6|) for $o(.^, U2(^)), we write the left hand side of the in- 
equality of (|3.2U|) as half of 

-V{0 + V{u2{t)) (3.21) 

-- y - - U2) [ , , , ^=]di ■ 

Ju2 (4 - /U)a/(Ui - /i)(/i - 1X2) 

The first two terms -V{0 + V{u2{t)) < -V{u2{t)/2) + V{u2{t) < for -1 < 
^ < U2{t)/2. The whole third term, in view of the scaling ()3.14|) . has an order in t 
lower than that of —V{u2(t)/2) + V{u2{t)). These prove that ()3.21|) is negative for 
— 1 < ,^ < M2(t)/2 when t is negatively large. We have therefore proved (j3.2(Jj) . 

Inequalities (j2.3UH2.3'T|) have thus been verified. By Theorem 2.3, the equilibrium 
measure for the external field V{^) = t^" + K(0 has no gap in its support when t is 
negatively large enough. 

The scaling nature of the above approach allows us to extend the result from 
p(0 — to an arbitrary monic polynomial of degree n. 

Theorem 3.1. Under conditions \3.8\3.lO^) on K(0> equilibrium measure for the 
external field V{^,t) = K(0 + ^ PiO> where p{^) is a monic polynomial of an odd 
degree n, has no gap in its support when t is negatively large enough. 
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The case when t is positively large can be treated in the same way. Here, we 
replace conditions ()3.8II3.11|) on K(0 by analogous ones 

K(0 = C.\^\''- + h40 fore«-l, (3.22) 

K(0 > fore»l, (3.23) 

where C_ > and M_ > max{2,n}. Function /i+(0 has the following behavior 



lim = . (3.24) 



Theorem 3.2. Under conditions II S. 2S9^3~24\ ) on K(0> equilibrium measure for 
the external field V{^, t) = K(0 + ^ P{0> where p(^) is a monic polynomial of an odd 
degree n, has no gap in its support when t is positively large enough. 

§ 3.2 The degree of polynomial p(^) is even 

We shall show that the equilibrium measure for the K(0 +^ P(0 has a single gap 
in its support when t is negatively large. When t is positively large, there are many 
possibilities. The equilibrium measure has no gap when p{S,) is convex and may have 
multi-gaps when p(^) is non-convex. 

We will first study the case when t is negatively large. The minimizer will be 
supported on [^4,^3] U [m2,Mi]. This corresponds to (7 = 1 in ()2.5|) . 

For simplicity, we assume that K(0 of (j3.1|) is an even function. The evenness 
of V{^) = t^" + K<(0 implies the evenness of the equilibrium measure V^(0- Conse- 
quently, its support must be symmetric about the origin. This means that U3 = —U2 
and M4 = —Ui. We therefore obtain from ()2.28|1 that 

P(e,^) = 2{e~u\){e-ul)^^{^,u) + {^ + Ur){e-ul)^i{ur,u) 

+ - + U2)<^1{U2,U) + (^2 _ - M2)^l(-U2, U) 
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Equations ()2.27|) exactly become 

2ui{ul — ul)'^i{ui,u) ui = 

71 

2U2{ul - ul)'^l{u2,u) U2 = 

which are equivalent to 



{u,-U2)^/2{u,+U2)^ ^^'^'^'^'''^ - = 0, (3.25) 

+ ^i(u2,m) = 0. (3.26) 

We have used the following identity when deriving ()3.25|) 

d 

^i(mi,m) - ^i(m2,m) = 2(mi - U2)^—'^i{ui,u) . 

OU2 

Its proof is the same as the one for a similar identity (j3.6j) . 

In view of the integral formula ()B.10|) for its left hand side, equation (j3.26j) is 
equivalent to 

^'(^'^^ -0. (3,27) 



U2 



We further assume K satisfies the condition 

K(0 = C|^l'' + MO for 1^1 »1, (3.28) 

where C > and M > max{3,n}. Here h{C,) satisfies 



Tm^ = ■ (3.29) 



This implies 



hm ^ = , hm ^ = , hm ^ = , lim = . (3.30) 

We will solve equations ()3.25ll3.25j) for large ui and U2. 

Again, we will first consider p{C,) = 

We split \E'i of (j3.25ll3.25|) into simpler terms 
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where and \l/gA/ are homogeneous functions of ui, U2) of orders n — 2 and M — 2, 
respectively. 
Introducing 

t/, = ^, £/. = -53, (3.31) 
we scale equations (|3. 25113 .251) as 

(t/i - t/2)V2(t/i + f/2){-^^e(f/i, [/i, ^72, -t/2, -t/i) 

^1' ^2, -f/i) (3.32) 
H m33-^*/i(0(^1'^1'"2, -M2, M = > 

|^|M-nCU2 ^|t|2(M-n) 



H 4^2:^/1(0(^1' "1' ^2, -M2, -Ml) - ^e(f^2, f/l, t^2, -t^2, "t^l) (3.33) 

\t\ ^'--^ 

+C^^M(f/2, f/l, f/2, -t/2, -f/l) H 432:^M0("2,Ml,M2, -M2, "Ml) = . 



At t = —00, equations ()3. 32113. 33|) become 

(f/l - f/2)V2(f/i + t/2){--^^e(t/l, f/l, f/2, -t/2, -f/l) 
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+C— -%.([/i,f/i,f/2,-f/2,-f/i)}5 = , (3.34) 
0U2 

-^e(f^i' f^2, -t/2, -f/l) + C7^^M(f/i, f/l, f/2, -f/2, -f/l) 

-^e(f^2, f/l, f/2, -f/2, -f/l) + C%.(f/2, f/l, f/2, -f/2, -f/l) = . (3.35) 
They have a solution Ui = U2 = U , where If is defined by 

-^^4U, U, U, -tj, -U) + C%.(f/, f/, U, -tj, -U) = i^K" + ct%=u = ' 

where we have used formula ()B.9|) for \I^^n and \E'gA/ in the first equality. We hence 
obtain 

as a solution of 1)3.34113.3,^ . 



It is straight forward to use the Imphcit Function Theory to determine the solution 
of equations ()3. 32113. 3^ in the neighborhood of Ui = U2 = U , t = —00. This in turn 
gives solution ui{t), U2(t) of ()3. 2 5113 .2^1) for t ^ — 1. The solution has the asymptotics 

ui{t) ^ U\t\M^ ^ ui{t) ^ U\t\M^ for t < -1 . 

We now verify inequalities (12.30112.311) . We will first prove ()2.3()j) and the last 
two inequalities of ()2.3H) . It suffices to show that ^i{C,,ui(t),U2{t), —U2it), —ui(t)) 
is negative for ^ < —U2{t)/2 and positive for ^ > U2{t)/2 when t ^ —1. Since 
'Ui(t), 'U2(t), — 'U2(t), — 'Ui(t)) is odd in ^ in view of formula ()B.8|1 . we only need 
to prove its positivity for ^ > U2{t)/2. 

We split $1 into simpler terms and scale these terms 

$l(^,Ml,M2, -U2, -Ml) 

= f/i, U2, -U2, -Ui) + f/i, U2, -U2, -Ui) 

+ ^Zi^h{i){.i,Ui,U2, -U2, -Ml)] , (3.36) 

f M-n 

where S = ^/|t| M-n . in view of condition ()3.3()|1 and the scaling 1)3.3111 . the last term 
in the parenthesis goes to zero uniformly for U/3 < C,/\t\ < 3U as t ^ — 00. The 
sum in the parenthesis then has the limit 

-$5^(2, U, U, -U, -U) + C$5m(H, U, U, -U, -U) 

2f/(S2 - f/2) 

where we have used ()B.7|) in the first equality and U = [n/CMY^^^~^ in the last one. 
The limit is positive for t//3 < H < 3U. It then follows from equation ()3.36p that 
^i{^,ui{t),U2{t),-U2{t),-ui{t)) is positive for U2{t)/2 < ^ < 2ui{t) when t < -1. 
We now show that $i(^, 'Ui(t), 'U2(t), —U2{t), —Ui{t)) > uniformly for ^ > 2ui{t). 
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We use formula ()B.8|) to write $i(^,mi,M2, —U2, —ui) as 



, 1 r Vif^'-^iw-^^i) v'if^)df^ 





(3.37) 




where we have used (|3.27|) in the equahty. 

It is then straight forward to use the argument below (j3.19p to show that the 
integral of ()3.37|) tends to zero uniformly for C, > 2ui(t) as t ^ — oo. 

The first term V'{^) = tn^^~^ + ^^e parenthesis is bounded from below 

by a positive constant uniformly for ^ > 2ui{t) when t <^ —1. This follows from the 
conditions ()3.28ll3.30j) on 

Function ^i{C,,ui(t),U2{t), —U2(t), —Ui(t)) is therefore positive in view of ()3.37|) . 
This together with the similar result for U2{t)/2 < C, < 2ui{t) proves that it is positive 
uniformly for ^ > U2(t)/2 when t ^ — 1. 

We have therefore proved ()2.30|) and the last two inequalities of ()2.31|) . 

It remains to prove the rest of ()2.3ip . i.e., 



for 1^1 < U2(t) when t < -1. 

First, the inequality ()3.38|) is vahd for —U2(t) < C, < —U2(t)/2 since $i is negative 
for —2ui{t) < C, < —U2{t)/2. It suffices to prove the inequality for —U2{t)/2 < ^ < 

U2{t). 

We use formula ()B.8|) to write the integral of ()3.38p as 

^(0 - V{-U2) 




V(/X2-Ml(t)')(/X'-M2(t)2) M2(t), -U2{t), ~uiit))dix > (3.38) 



2 



1 

+- 



/ sj{u^-ul){u^-ul)[l 



V\fi)dfi 




TT 
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It is easy to use the scaling ()3.31|) to show that the above is positive for —U2(t)/2 < ^ < 
0; hence, inequahty ()3.38|) is verified for —U2(t) < ^ < 0. Since -ui, ^2, —^2, — mi) 
is odd in in view of ()B.8j) . the inequahty ()3.38|1 can be extended from —U2(t) < ^ < 
to —U2{t) < ^ < U2{t). This completes the verification of inequalities (12.30112.311) . 

By Theorem 2.3, the equilibrium measure for the external field V{E,) = tC" + K(0 
has a single gap in its support when t is negatively large. 

It is also easy to extend the result from p(^) = ^" to any even monic polynomial 
of degree n. 

Theorem 3.3. Suppose K(0 ^■^ on even function that satisfies condition \3. i^<jCO^) . 

The equilibrium measure for the external field V{C,,t) = K(0 + PiO' where p{^) 
is an even monic polynomial of degree n, has a single gap in its support when t is 
negatively large. 

For the case that t is positively large, we assume 

K"(0 > for lel » 1 . (3.39) 

Again, we first consider the case p(^) = 

Since is uniformly bounded in |^| < A for some large constant A, we may 

choose t positively large enough so that V{^) = t^" + K(0 is a convex function for 
1^1 > Co 5 where eo is a tiny positive number. It is well known that an everywhere 
convex external field has an equilibrium measure whose support is a connected finite 
interval [20j. Hence, if V{^) = t^" + K(0 is also convex in |^| < eo for large t > 0, 
the corresponding equilibrium measure has no gap in its support. 

It is therefore interesting to note that, even if V{C,) = t^" + K(0 is never convex 
in a small neighborhood of ^ = 0, the equilibrium measure can still be shown to have 
no gap in its support for large t > 0. 

The approach is similar to one we use in § 3.1. We will present the proof briefiy. 

We still need to solve equations ()3.2H3.3|) . Using 

Ui = Uit^ , U2 = U2t" , 
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we scale the equations as 



{Ui - U2)J^^^r.{U,,U^,U2) + ^j^^v,iuuu^,u2) -J- = 0, (3.40) 

V OU2 OU2 V VT 



^^n{U,,U,,U2) 

1 1 

+ — ^K(Ul,Ui,M2) + ^e(^2,^7l,f/2) + — ^K(?i2,Ml,M2) =0 . (3.41) 

At t = +00, these equations become 



(f/l-f/2)'^^%>(f/l,f/l,t/2)-y^ = 

^^4U^,U,,U2) + ^^n{U2,U,,U2) = 
They have a solution Ui = f/, f/2 = — f/, where U is defined by 



It is not hard to get from a formula of type ()B.11|) for \E'gn that 



dU2 ' ' ' ' ' 4(n-2)!! 

We hence obtain 




as a solution of equations 

One can then use the Implicit Function Theory to show that (|3. 40113. 41}) give IJ\ 
and f/2 as functions of t near t = +00 and 

limt^+ooUi{t) = U , limt^+ooU2{t) = ~U . 

Therefore, equations ()3.2ll3.3p can be inverted to give ui(t) and U2it), which have the 
asymptotics 

Ui[t) ^ — , U2[t) ^ — Y as t > 1 . 
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To verify the inequahties ()2.30II2.3T|) . it suffices to prove that $o(^, ui(t) , U2{t)) > 
uniformly for all ^ when t is sufficiently large. 
We first decompose $0 

= tt[$^„(4,f/l(t),f/2(t)) + 4'^'K{C)(^>«lW>«2(t))] . 

tn f,n 

The first term in the parenthesis is bigger than a positive constant for all ^ and 
large t. To see this, function $gn, in view of formula ()B.4|) . can be written as a 
multiple integral of the second derivative of S,^. Hence, when U/2 < Ui < 2U and 
—2U < U2 < —U/2, the first term is bounded from below by a positive constant for 
all e/t^/". 

To show $o(^, ui(t) , U2{t)) > for all ^, it is enough to show that the second term 
^v*{^) is bounded from below for all ^ and large t. To accomplish this, we deduct 
from condition (j3.39|) that V"{C,) is bounded from below for all ^. Function $v,(5); 
which can be written as a multiple integral of V" in view of formula ()B.4|) . is therefore 
bounded from below for all ^, Ui and U2. 

We have therefore verified inequalities (j2. 30112 .311) . By Theorem 2.3, the equilib- 
rium measure for the external field V{C,,t) = K(0 + '^^^ is supported on a single 
interval when t ^ 1. 

To generalize p{^) from ^" to any convex polynomial p{C), we notice that the 
lowest order term ^"^ in p(^) must be of an even order and that its coefficient must 
also be positive. The previous analysis centered around can be applied to 

Theorem 3.4. Under condition \3. 39^) . the equilibrium measure for the external field 
V{C,,t) = K(0 +^ P{0> where p{^) is a convex polynomial, has no gap in its support 
when t ^ 1. 

We conclude this section by making an observation on the case of non-convex p{C,)- 
Such a polynomial p(^) has multiple "wells". They will be amplified and become the 
"wells" of the external field V{C,,t) = + ^ PiO ^ is positively large if K(0 
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satisfies condition ()3.39p . Its equilibrium measure is therefore likely to be supported 
on multiple disjoint intervals. 

Appendix A. Algebro- Geometric Solution of the Riemann- 
Hilbert Problem 

In Section 2, we use function theoretical methods to solve the Riemann-Hilbert 
problem (j2.14j) and (j2.15|) . Our solution is given by formula (|2.16p : it is the cor- 
nerstone of Propositions 2.1 and 2.2. In this appendix, we will present yet another 
approach to the Riemann-Hilbert problem. We will solve it using an algebro-geometric 
method. More precisely, we will give another expression for the solution ()2.16|1 and 
hence re-derive the formulae of Propositions 2.1 and 2.2. 

The Riemann-Hilbert problem ()2.14|) and ()2.15p has an intrinsic algebro-geometric 
structure. All its solutions are connected to the Riemann surface defined by the 
equation w'^ = {fi — ui){fi — U2) ■ ■ ■ (/i — U2g+2)- We choose the branch cuts along 
the set of / of ()2.5|1 . It is remarkable that our algebro-geometric approach does not 
require the external field V{S,) to be an analytic function. 

Our starting point is to generalize the Cauchy kernel. On a Riemann surface, there 
are many analogues of the Cauchy kernel. The most convenient one for the Riemann- 
Hilbert problem ()2.14j) and ()2.15|) is an Abelian differential of the third kind, denoted 
by K{fi,z)d^, with two simple poles at the points {z,±R{z,u)) with residues ±1, 
respectively. Here, ±R{z, u) are the upper and lower sheets of the Riemann surface. 
This means that K{^^ z)dfi takes the form 

T^/ \ 7 R(z,u) 

A f/i, z)a^ = — — ^ h holomorphic terms 

it(/i, u) n — z 

and that its behavior, as ^ is near the poles, is 

K{^, z)dfj. = — [- regular terms. 
{X — z 

We may further require that 

K{ix,z)dfx = k = l,2,--- ,g , (A.l) 
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T^, X , '^/^ R(z,u) , , , r 

K{^,z)d^= — — — + > / 

R{fi,u) f^-z f-^ A, 



which is analogous to condition ()2.13p . 

The differential K{fi, z)dfi then takes the form 

dfi R{z,u) , , R{z,u) 
k=i ou,,+,R{i.u) i-z 

where ujk{fi)dfi is the basis of holomorphic differentials normalized along the intervals 
[u2k+i,U2k], k=l,...,g. 

Using the Riemann bilinear relations between K{fi, z)dfi and the differentials 



-dr] defined in fl2.22|) . it is possible to reduce the above formula to the form 



-R{f^,u)^^^ 



where p± = {z,±R{z,u)) and Tk{u) are the coefficients of the expansion (j2.25p . 
Formula ()A.2|) can also be derived from the explicit form of the Bergmann kernel on 
the Riemann surface w"^ = R^{^,u) 

We next point out a remarkable symmetry property 

d d 
-Ki^,z) = -Kiz,f.). 

It also follows from the Riemann bilinear relations 

We now use K{fi, z)d^ to construct a solution of the Riemann-Hilbert problem 
and (ITTH|1 

g{z) = -- [ V'{^,)K{z, + ^P^^o(z,u) 

IT J I IT R{Z,U) 

= TT ~ / V{f^)K{fi,z)dfi + p. • 

OZ ITT J J J ^ R{Z, U) 

To see this, we derive from ()A.3|) the boundary value of Q at real ^ 



(A.4) 



G+{i) = v^y'(0 - -P.V. I V'{fi)K{^,fi)dfi + 

Jl 



When ^ G / and fi & I, the kernel K{^, /i) is real and R{C,, u) is pure imaginary. We 
then derive ImQ^{^) = V'{^), which is fl2.14p . When ^ G M\/ and /i G /, the kernel 
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K{^,fi) is pure imaginary and R{^,u) is real. We instead have ReQ^{C,) = 0, which 
is exactly ()2.15|) . 

We further claim that Q of ()A.3j) equals the function defined in ()2.1(j|l . To see 
this, it suffices to prove that Q of ()A.3|) satisfies both the loop conditions ()2.13|) and 
the asymptotics Q{z) = — +0(1/2;^) for large \z\. Conditions ()2.13|) are easily 
verified, in view of similar conditions ()A.1|) on K{fi,z)dfi and ()2.23|) on PgQ. For 
large \z\, the kernel K{z,fi) behaves like O(^). Indeed the first term of K{z,fi) is 
proportional to u){z-n) ' which clearly decays at least as fast as as \z\ — > oo. 
The other terms of K(z^fi) are of the form ^^^^ = 0( ^J+i ) for \z\ large and m = 
l,...,g. This, together with the behavior (j2.22j) of Pg^, justifies the asymptotics 
g{z) = + 0(1/^2) for large \z\. 

Inserting the explicit form ()A.2|) into the expression ()A.4|) . we obtain 



R{z,u) f V{n)djj, 



IT J J R{jj,, u){ijl — z) 



(A.5) 



+ 2V ^2^^Qa,m{u)2_^^kL^_,[u) + ^ , 

where qg^miu) is given in ()2.2(')|1 . 

The boundary value of Q {z) of ()A.5|) on the real axis can be obtained by observing 
that the first term has the boundary value at real ^ 

^ [-2v^i?(e, u)mg{i, u) + v^y(0] , 

where is as given in Proposition 2.2. The function '^g{^,u) is related to (^g{z,u) 
of Proposition 2.1 by an identity J2] 

(A,6) 

We hence arrive at the boundary value of Q at the real ^ 

g+^^. ^ -2V^R\i,u)^g{i,u) + ^Q{i) ^ yzi^'(^) 

R{i, u) 
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where the polynomial Q{^,u) is 

29+2 ^ , g m ^ 

i=l ^ * m=l k=l 

These are equivalent to the formulae of Propositions 2.1 and 2.2. 

Appendix B. Euler-Poisson-Darboux Equations 

The boundary value problem for the Euler-Poisson-Darboux equations 
has one and only one solution. Its solution can be constructed using those of the 
following simpler problem as building blocks [211 1121 



2{xi - X2 



dxidx2 dxi 3x2 
g(xi,xi) = g{xi) . 



dq dq 

p— — , p > ts a constant 



(B.l) 
(B.2) 



A simple calculation shows that the solution of ()B.HIB.2j) is given by the formula 



q{xi, X2) = Co / , (1 + At) 2 dp , 



p- 



where 



Co 



■1-1 ^/iZJp 



A change of integration variable gives another formula for the solution 
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q{xi,X2) = Co 



p-i 

~2~ r^i 



AX-X2) 2 

TT^^ 

X2 [Xi — X) 2 



_Xi - X2_ 

where the square root is set to be positive for x between xi and X2- 

Using a solution method of [21], one is able to construct the solution of 
using a multiple integral|jl2j 

-1 



= Mo 



2 



/ ^ ^ ^ ^ - dpi ■ ■ ■ dp2g+2 



1 



nr=f(i-/^.) 



where / denotes 
/ 



(g+2) A +/^2g+2 



(B.3) 



(B.4) 



.1+ P2A + ^il I- Pl . . 1 - ^i2g+2 . 

(^— ^ + ^—^1) + ■■■) + — ^^^23+2) 
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and the constant Mq is chosen so that the boundary condition ()2.20|) is satisfied. Here 
j(9+2) (denotes the {g + 2)th derivative of /. 

The function \t'c,(^,?2) satisfies the same equations ()2. 18112. as ^g{C,,u); but, 
instead of the boundary condition ()2.20|) . it satisfies ()2.24|) with one derivative lower. 

We now hst some of the properties concerning $0) "^i s-nd They are useful in 
the calculations in Section 3. 



Property B.l. 



V(^-^l)(^-^2) 

1 n V'{fx)dfx 



2^ Ju2 - u) a/(ui - u){u - U2) 



(B.6) 



2(4^ — -U^jM 



^l(^,Ml,M2, -U2, -Ul] 



2V(e-«?)(e-«i) 



^Ju2 {^i^-e)^/W^W)W^ 
'$i(u,u,u,-u,-u) = ^ , (B.9) 

1 r(/i)c;/i 



, = . (B.IO) 

Function V{^) is assumed to be an even function in ^B. 7| ), iB.8\} . HB.9\} and HB.10\} . 

Proof. In formula ()B.4|1 . $0) '^'i and \E'i are written as multiple integrals of V"{C,), 
V"'{C,) and V"{^), respectively. Since the smooth function V"(^) can be approximated 
by polynomials in C^{S) on every compact set S as close as possible, it suffices to 
prove Property B.l when V{^) is simply a polynomial. 
We will rely on the contour integral formulation of 

t,K.s) = -^/^i»-. (B.11) 
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where 7 is a contour enclosing the point ^ and all the cuts along [u2k,U2k-i], k = 
1,2, - ■ ■ ,g + l. To see this , it is easy to check that the right hand side, as a function 
of ^ and u, satisfies fl2.18j) and ()2.19|1 . Boundary condition ()2.2()j) is also easily verified 
using the Cauchy Integral Formula. 

Identity ()B.5|) is an easy consequence of ()B.11|) when ui = U2 = u. 

To prove ()B.6j) . we replace the contour 7 by 7', which still encloses the cuts [u2, ui], 
k = 1,2, ■ ■ ■ , g + 1, but excludes the point ^, and write $0(^5 wi, ^2) as 



The first term gives the first term on the right of ()B.6|) . Rewriting the second term 
as an integral along the cut [^2, Ui] gives the second term of ()B.6|) . 
Identities ()B.7|) and ()B.8|) can be proved in the same way. 

To prove ()B.9|1 . we use an analogous formula of ()B.11|) for Ui, U2, —U2, —Ui) 



Identity ()B.9|) is an immediate consequence of this formula and the Cauchy Integral 



To prove ()B.10|) . we notice from the formula of type ()B.4|) that \t'i can be writ- 
ten as a multiple integral of V"{E,)- Since V{S,) is an even function of ^, so is 
'^i{^,Ui,U2, —U2, ~ui). This evenness allows us to modify formula ()B.12|) to get a 
new contour integral formulation for \E'(^, Ui,U2, —U2, —Ui) 




(yU - y/{fI-Ui){n-U2)\ 47r Jj' (/^ - 0^(/^, ^) 





(B.12) 



Formula. 




The left hand side of flB.10|) then equals 
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We obtain the right hand side of ()B.10|) by writing the last integral along the cuts 
[—Ml, — M2] and ['U2,'Ui] and using the oddness of V'{fj,). 
The proof of Property B.l is completed. 
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